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Description 

METHOD FOR ADAPTIVE SEGMENT 
REFINEMENT IN OPTICAL PROXIMITY 

CORRECTION 

Background of Invention 

[0001] The present invention relates in general to lithographic 

processes and, more particularly, to a method for correct- 
ing a mask layout of an integrated circuit for lithographic 
processing effects. 

[0002] The advent of advanced lithographic techniques and the 
attempt of chip manufacturers to closely follow Moore's 
law predicting an exponential growth of number of com- 
ponents on a chip and their shrinkage on the wafer is 
making the process of designing tools for automating the 
chip design extremely challenging. The chip is imprinted 
by means of optical lithographic techniques on the silicon 
wafer by way of chrome on glassmasks. As the compo- 
nents on the chip become smaller, they are now in the 
verge of reaching the limits of capacity of the lithographic 



process. 

[0003] one of the most common problem of printing integrated 
circuits on silicon wafers from the mask using lithographic 
processing is the loss of fidelity of the mask shapes. The 
dimensions of the wafer shapes are usually much smaller 
than the wavelength of light that is used in the litho- 
graphic processing. The non-linearities associated with 
this lithographic process of imprinting and the laws of 
physics associated with light (including diffraction of light 
waves) makes impressing highly unpredictable. For in- 
stance, rectangles are foreshortened into elliptical- 
shapes, right angles are rounded, and the width of the 
shapes shrinks. In addition to the non-linearity associated 
with the optical effects, the process of developing after 
exposing the wafer and the photoactive material (i.e. pho- 
toresist, equivalently referred to as resist) thereon also 
adds to the uncertainty associated with the critical dimen- 
sions of the features. Other effects mainly related to the 
exposure and development effect associated with the re- 
sist exist that contribute to the distortions. The wafer 
shapes get deformed due the diffraction and scattering of 
light and other related optical and chemical effects of 
photo-resist materials. The combination of these litho- 



graphic processing effects are known as optical proximity 
effects. 

[0004] jo counter the problem of optical proximity, and increase 
the fidelity of the wafer printing, mask designers take into 
account the deformities introduced by the aforementioned 
lithographic processing effects by intentionally and sys- 
tematically distorting the original shapes on the mask. 
The net result of these distortions is that the imprinted 
shape on the wafer ultimately looks like the target or in- 
tended images, satisfying the design rules that were cre- 
ated to achieve the desired yield in chip manufacturing. 
These methods are generally referred to as optical prox- 
imity corrections (OPC) and can be categorized into three 
classes: 

[0005] Ad-hoc Method: This method is almost as old as manufac- 
turing VLSI chips. Early designers modified existing meth- 
ods by putting "flares" and "hammer heads"at the end of 
lines and "serifs"at the rectangular corners to compensate 
for the "line-end shortening" and "corner-rounding" due 
to process irregularities. Since at the early stages in the 
history of chip manufacturing, the size of the features was 
still large compared to the wavelength of the light used in 
the lithographic process, optical proximity effects did not 



contribute significantly to the total error budget of chip 
manufacturing. However, as the size of the features have 
continued to shrink, mask designers have continuously 
readapted earlier techniques to obtain the desired results. 
Therefore, though prevalent in the early part of the history 
of optical proximity effect corrections, there has not been 
any significant development of these so-called ad-hoc 
methods. 

[0006] Rules-Based Method: The rules-based method is an initial at- 
tempt to formalize the above ad-hoc method. It was ob- 
served by chip designers that any compensation required 
by a particular shape on the wafer was dependent on 
neighboring shapes, e.g., an isolated line requires more 
compensation than a set of dense lines. The rules-based 
method formalizes this notion in a more quantitative way. 
For example, a chip designer may use a set of empirically- 
derived rules to decide what the dimensions of a "ham- 
mer- head" are to be added to compensate for line-end 
shortening or what "corner-serifs" would compensate for 
corner-rounding. These dimensions can be applied as a 
function of certain rules as, for instance, those that de- 
pend on the particular shape to which it is applied, the 
distance to the neighbors, and the dimensions of the 



neighbors. However, the rules are still applied heuristi- 
cally, and can still lead to inaccurate imaging relative to 
the desired target image, even when implemented in a 
computer program. 
[0007] Model-Based Method: This method, also known as model- 
based OPC (MBOPC), emulates the physical and optical ef- 
fects that are mostly responsible for shape deformations. 
At the heart of these methods is a computer simulation 
program that, given the appropriate optical and physical 
parameters and the original dimension of the object on 
the mask, predicts, with a certain degree of accuracy, the 
printed dimension of the object on the wafer. In the cor- 
rection phase of the model-based optical proximity cor- 
rection (MBOPC), the shape on the mask is iteratively 
modified so that the resulting output closely approxi- 
mates what is desired for the image or imprinted shape 
on the wafer. This method automatically deforms existing 
mask shapes to achieve the target dimensions on the 
wafer. 

[0008] The success of model-based OPC depends on a very ac- 
curate simulator that would predict the lithographic pro- 
cessing effects. The simulator predicts the lithographic 
effect corresponding to selected points on the mask. 



Based on the simulation, an integrated circuit layout mod- 
eler can determine the printed image on the wafer by in- 
terpolating selected simulated points. The modeler may 
also correct the mask shape by modifying the edge seg- 
ments on the mask based on the simulation done at rep- 
resentative points of a mask segment, typically on the 
center point of the edge segment. The simulation of the 
image at points on the mask is computationally very ex- 
pensive. Therefore, the modeler usually relies on doing 
the simulation only at certain predetermined or prese- 
lected evaluation points. The preselection is done by seg- 
menting the edges of the mask shapes by using certain 
pre-specified rules. 
[0009] The preselection of points often make the nature of simu- 
lation as a whole unreliable and untrustworthy. This is be- 
cause often the detailed variations of the wafer image 
cannot be characterized by simple interpolation of these 
simulated images done at the preselected evaluation 
points. The simulation at the preselected segment points 
often has the problem that more points of simulation than 
necessary are provided where there is little image varia- 
tion at the wafer, while having too few points of simula- 
tion where there are large image variations. 



[0010] The current state of the art is shown in Figure 1. In the 
current state of the art, the input is a mask layout (101) 
that is intended to project a desired target image 106. An 
example of an input mask layout 101 is shown in Figure 
1A. The mask shapes are segmented at a preprocessing 
stage (Block 102), resulting in segmented mask shapes 
103. The resulting mask segmented mask shapes 103 is 
similar to that illustrated in Figure IB. This stage precedes 
any image simulation (Block 104). After the edges are 
segmented, the image intensity is simulated typically at a 
predefined point within each segment only (Block 104). 
The image at that predefined point of the segment is then 
compared (Block 105) with the corresponding points of 
the target image 106. (An example of a target image 106 
is shown in Figure 1C.) If the difference of image intensity 
at the predefined point of the segments and the corre- 
sponding points of the target image 106 is greater than 
the specified tolerance, segments of the mask are moved 
away or towards the target shape 106 to compensate for 
the difference (Block 107). This process is repeated until 
all the intensity values at the predefined evaluation points 
of all the segments are within the specified tolerance. At 
the end the output mask layout 108 is produced. For clar- 



ity, the initial mask shapes 101 are illustrated to be the 
same as the target shapes 106. Initially, the mask shapes 
101 may be the same as the target shapes 106. However, 
after a few OPC iterations (e.g., Blocks 104, 105 and 107), 
the output mask shapes 108 become different from the 
target shapes 106. An example of the output mask shape 
108 is illustrated in Figure ID. 
[001 1] jo accelerate the process of convergence during the itera- 
tive optical correction process, as described above, the 
amount of movement of the mask segments at each itera- 
tion is controlled by the gradient and curvature of the im- 
age intensity profile at the vicinity of the above mentioned 
predefined evaluation point for each segment. For a given 
value of intensity l(x, y) at a point (x, y) on the wafer, the 
gradient and the curvature is defined as follows: 



[0012] where l(x, y) is the image intensity at the predefined point, 
and where, in this example, the x-direction extends along 
the directions of segments A and B, and the y-direction is 
the direction orthogonal to the segments. In order to 
compute the gradient accurately, a finite difference 



The gradient is defined as 




(1) 



method is used in the current state of the art tools. In this 
method, in addition to the computation of the image in- 
tensity at the selected point, image intensities at a few 
more points are also computed. An example of this ap- 
proach is shown in Figure IE. In Figure IE, a segmented 
layout is shown (103). One of the segments (153) has 
end-points 151 and 152, respectively. The image inten- 
sity, as explained above is computed at a predefined eval- 
uation point 160 (preferably the center) for segment 153. 
The image computed at the preselected evaluation point 
160 is assumed to be representative of the segment 153. 
In addition to the image intensity computed at the prede- 
termined evaluation point 160, the image intensities are 
computed at a few more points in the vicinity of the eval- 
uation point 160. Referring to Figure IE, the region of 
vicinity is marked by a square box 155 centered around 
the evaluation point 160 with edges parallel to the x and 
the y axis (note that for the sake simplicity, it is assumed 
that the edge 153 is parallel to the x axis). The size of the 
box 155 is very small compared to the length of the seg- 
ment 153. For example, if the length 171 of the segment 
153 is about 100 nanometers, then the length 170 of a 
side of the square box 155 (i.e. the vicinity of the evalua- 



tion point 160) can be as small as 2 nanometers. The im- 
age intensity at the given (predefined) evaluation point 
160 is defined as 10. The image intensities of two more 
points one immediately above and another immediately 
below (points 161 and 162, respectively) are defined as il 
and 12, respectively. Similarly, the image intensities of two 
more points one immediately to the right and another im- 
mediately to the left of 160 (points 163 and 164, respec- 
tively) are defined as 13 and 14, respectively. The distance 
between the points above and below the given point 
(along the y-axis) are defined as Ay , and the distance be- 
tween the points left and right of the given point (along 
the x-axis) is defined as Ax . Then the gradient is approx- 
imated as: 



The curvature is defined as: 
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[0013] The curvature can also be computed by a similar finite 
difference computation with a few more extra points, as 
would be understood by one skilled in the art. 

[0014] | n order to compute the curvature one need to consider 

points 165, 166,167,168 around the point 160 in the box 
155, along with the points 161,162,163,164 as previously 
considered for computing the gradient. The image inten- 
sities at these four points are computed as 15, 16, 17, 18, 
respectively. In this case the curvature may be computed 
as: 



[0015] Note that the above expression in Equation (4) is only one 
way of approximating the curvature by a finite difference 
method. Other methods may be used to approximate the 
curvature. For example, a simpler but less accurate 
method of computation of gradient using only the afore- 
mentioned five points 160 through 164 would be as given 
in the following equation 




(4) 



(5) 



[0016] | n a practical simulator the value of Ax and Ay is on the 
order of a few nanometers. It should be noted that the 
values of gradient and curvature are computed in current 
state of the art methods and are thus available for use in 
the context of the present inventive method and will not 
require any additional computation for the gradient and 
curvature to implement the present inventive method. 

[0017] T ne problem with the current method is demonstrated in 
Figures 2 through 5. Figure 2A shows a target image 200 
and an input mask shape 201 (which initially has the same 
shape as the target image 200) and the corresponding 
output mask shape 208. Note that the target image shape 
boundary 200 represents the contour of threshold image 
intensity that is desired in order to print the image. Here- 
inafter, references to the boundary and shape of an image 
are used to interchangeably refer to the contour of the 
threshold image intensity. It also shows the corresponding 
simulated image 220 (that is, the contour of the aerial im- 
age at the threshold intensity required to print the image) 
on the wafer and then finally superimposes the wafer im- 



age 220 on the desired target image 200. The image in- 
tensity threshold level is defined by the corresponding 
lithographic process conditions, including, but not limited 
to, for example, factors such as the dose-to-clear charac- 
teristics of the resist. Figure 2 shows a portion of a target 
shape 200 having (straight-edged) boundary 215, which 
has been divided into edge segments A (having segment 
boundaries 210, 212), B (having segment boundaries 212, 
214), etc. Predetermined evaluation points (e.g. points 
310, 320 in segments A, B respectively) are defined for 
each segment, which are typically chosen to be at the 
center of a segment. Superimposed on the target shape 
200 is the corresponding (simulated) wafer image 220 
having (curved-edged) boundary 225. The interior of 
boundary 225 of the simulated image intensity is defined 
as the region where the image intensity exceeds a partic- 
ular threshold intensity level. The current state of the art, 
the image intensity is typically computed only at the cen- 
ter or some predefined point of the segments A, B, etc. As 
described above the image gradient and the curvature is 
also computed along with the image intensity. The OPC 
algorithm incrementally corrects (Block 107) the mask 
shapes until the simulated image 220 has a threshold in- 



tensity that matches the target image intensity at the 
points of the target image corresponding to the predeter- 
mined evaluation or comparison points 310, 320 of each 
of the mask segments (See Figure 3). The conventional 
OPC correction algorithm ignores image intensities at any 
point that is between two pre-selected comparison points 
(e.g. points 310, 320) of adjacent target image segments. 
As a result of this, one can get severe deviations of the 
wafer image from the target image. An example of this is 
shown in Figure 4. 
[0018] FIG. 4 shows a target feature 400 having a straight-edged 
threshold boundary 410 that outlines a large shape 401 
adjacent to a narrow shape 403. The simulated image 411 
having a (curve-edged) threshold boundary 413 resulting 
from a mask generated by conventional MBOPC is super- 
imposed on the target image boundary 410. Because of 
the failure of the current state of the art it creates a 
neckingor constriction, wherein the width of the resulting 
image W at the neck region 402 is less than the critical 
dimension (CD) of the narrow target image shape 403. 
Sometimes the constriction can become really narrow and 
pinch-off as shown at the center of the scanning electron 
microscope (SEM) image in Figure 5. In FIG. 5, a wide fea- 



ture 501 (analogous to the wide portion 401 of the target 
image in FIG. 4) and a narrow feature 503 (analogous to 
the narrow portion 403 of the target image in FIG. 4) 
should have been connected at point 502, but in fact has 
printed as two disconnected shapes due to excessive 
necking of the image. In the conventional MBOPC, the so- 
lution to this problem requires creation of smaller seg- 
ments throughout the mask and target image, which will 
improve the accuracy of the simulation and allow the OPC 
model to modify the mask in the critical neck region. 
However, in conventional methods, the refined segmenta- 
tion is performed even in portions of the mask where the 
resulting image does not vary rapidly (e.g. along the nar- 
row region 213 away from the neck region 222 in Figure 
2A), which can significantly increase computation costs 
unnecessarily. 

[0019] Therefore, there is a need for a methodology for design- 
ing masks that result in images that match target shapes 

accurately, but without undue computation time. 
Summary of Invention 

[0020] it is therefore an objective of the present invention to pro- 
vide a method and a computer program product to per- 
form the method, for manufacturing an integrated circuit, 



including designing a lithographic mask layout that re- 
sults in a printed image that accurately represents a target 
image without significant increases in computation costs 
over conventional model-based optical proximity correc- 
tion algorithms (MBOPC). 

[0021] | t j S a further objective of the present invention to provide 
a method and computer program product for performing 
MBOPC in which mask segments are adaptively refined. 

[0022] | t j S another objective of the present invention to provide 
a method and computer program product for adaptively 
refining mask segments without significantly increasing 
the computation costs by basing the refinement on image 
computations at pre-existing segment evaluation points. 

[0023] a methodology is presented for a integrated circuit layout 
design tool to improve the accuracy and correctness of 
optical proximity correction. The methodology adaptively 
refines the segmentation of the mask shapes to improve 
the accuracy of the simulation. The adaptive refinement 
approach uses an iterative algorithm. Each step of the al- 
gorithm uses the local image gradient and image curva- 
ture to determine the deviation of the image from that 
predicted by the existing fragmentation. Instead of using 
a simulation of the exact contour it uses a curve fitting of 



the contour by using one or more of the parameters of the 
image intensity value, gradient and the curvature. If the 
deviation is more than a tolerable limit, it refines the ex- 
isting segmentation by inserting more segments at de- 
sired locations. The methodology can be used for im- 
proved simulation of the silicon wafer containing an inte- 
grated circuit and correct the mask that is used to create 
the silicon wafer. 
[0024] | n a fj rst aspect of the present invention, a method of 

manufacturing an integrated circuit, and a computer pro- 
gram product for performing the method, is provided in 
which a target design is provided, preferably expressed as 
a desired threshold intensity contour layout, wherein the 
target design shape is divided into target segments and 
each of the target segments has an associated segment 
evaluation point, providing an initial mask shape which is 
divided into mask segments corresponding to each of the 
target segments, determining a simulated image corre- 
sponding to each of the associated segment evaluation 
points, wherein values of the simulated image are formed 
in accordance with the mask shape for a given litho- 
graphic process, providing a deviation tolerance for a de- 
viation of a simulated image metric from a corresponding 



target image metric, and for each of the segment evalua- 
tion points, determining a simulated image metric based 
on the simulated image at said segment evaluation point, 
and refining the mask segment corresponding to the seg- 
ment evaluation point if the simulated image metric ex- 
ceeds the deviation tolerance. 

[0025] | n another aspect of the present invention, the deviation 
tolerance is the defined to be the maximum distance tol- 
erated between the threshold intensity contour of the 
simulated image from the target image. 

[0026] | n y e t another aspect of the present invention, an esti- 
mated image between segment evaluation points is used 
to determine the deviation metric. The estimated image 
may be obtained by a curve fit based on the simulated im- 
age at the evaluation points, such as intensity, gradient 
and curvature. A curve fit may be performed by any of a 
Binomial spline, polygonal approximation, circular arc, 
cubic spline, Bezier curve, or the like. 

[0027] in a further aspect of the present invention, pre-existing 
values of the gradient and curvature of the simulated im- 
age may be tested against a gradient and/or curvature 
tolerance to determine whether a segment refinement is 
required. If the gradient and/or curvature tolerances are 



not exceeded, then no further refinement of a segment is 
required, and no additional computation will be required. 
[0028] The foregoing has outlined rather broadly the features 

and technical advantages of the present invention in order 
that the detailed description of the invention that follows 
may be better understood. Additional features and advan- 
tages of the invention will be described hereinafter which 

form the subject of the claims of the invention. 
Brief Description of Drawings 

[0029] For a more complete understanding of the present inven- 
tion, and the advantages thereof, reference is now made 
to the following descriptions taken in conjunction with the 
accompanying drawings, in which: 

[0030] Figure 1 illustrates a flow chart for a conventional model- 
based optical proximity correction (MBOPC). 

[0031] Figure 1A illustrates an initial mask layout provided as in- 
put to MBOPC. 

[0032] Figure IB illustrates a mask shape that has been seg- 
mented according to conventional MBOPC. 

[0033] Figure 1C illustrates a target image pattern provided as 
input to MBOPC. 

[0034] Figure ID illustrates an output mask layout resulting from 
a conventional MBOPC. 



[0035] Figure IE illustrates grid points used in a finite difference 
calculation in a conventional MBOPC. 

[0036] Figure 2 illustrates a portion of a simulated wafer image 
superimposed on a corresponding portion of a target im- 
age having segments and predetermined evaluation 
points. 

[0037] Figure 2A illustrates an input target image and mask 
shape, a corresponding output mask shape, a corre- 
sponding simulated wafer image, and a superposition of 
the wafer image on the target image according to a con- 
ventional MBOPC. 

[0038] Figure 3 illustrates a portion of a simulated wafer image 
superimposed on a corresponding portion of a target im- 
age having segments and predetermined evaluation 
points, where the image intensity matches the target im- 
age intensity at the evaluation points. 

[0039] Figure 4 illustrates a portion of a target image feature 
having a non-critical shape adjacent to a narrow shape 
having a critical dimension (CD), and a superimposed cor- 
responding portion of a conventionally simulated wafer 
image superimposed on the target image, where the wafer 
image has a reduced width compared to the CD. 

[0040] Figure 5 shows an scanning electron microscope (SEM) 



image of a printed wafer image resulting from a mask lay- 
out output from a conventional MBOPC, where the wafer 
image has a pinched-off portion at the junction of a non- 
critical portion and critical portions of a target feature. 

[0041] Figure 6 illustrates a flow chart of a preferred embodi- 
ment of the present invention. 

[0042] Figure 7 illustrates a portion of a segmented target image 
having predetermined evaluation points, and a superim- 
posed corresponding portion of a simulated wafer image 
at a step of the method in accordance with the present in- 
vention. 

[0043] Figure 8 illustrates a portion of a segmented target image 
where a curve fit to the wafer image is performed between 
two predetermined evaluation points in accordance with 
the present invention. 

[0044] Figure 9 illustrates the deviation of the fitted curve from 
the target image in accordance with the present invention. 

[0045] Figures 10, 10A, 10B and 10C illustrate segment refine- 
ment in accordance with the present invention. 

[0046] Figure 11 illustrates a curve fitting performed on a refined 
subsegment in accordance with the present invention. 

[0047] Figure 12 illustrates a schematic diagram of a digital 

computer, adapted to perform the method of the present 



invention. 
Detailed Description 

[0048] | n the following description, numerous specific details 

may be set forth to provide a thorough understanding of 
the present invention. However, it will be obvious to those 
skilled in the art that the present invention may be prac- 
ticed without such specific details. 

[0049] Refer now to the drawings wherein depicted elements are 
not necessarily shown to scale and wherein like or similar 
elements are designated by the same reference numeral 
through the several views. 

[0050] | n accordance with the present invention, a method and 

computer system is provided for performing model-based 
optical proximity correction (MBOPC), wherein the com- 
parison of simulated image to target image is performed 
using an adaptive refinement of segmentation that is 
based on the image gradient. Using this methodology, a 
proper balance of the segmentation can be achieved that 
takes into account the variations of the wafer image. This 
methodology can be used to obtain more accurate and 
faster image simulation and mask correction tools. 

[0051] Figure 6 shows the preferred embodiment of the present 
invention. At the start, an initial mask layout is first pro- 



vided (Block 601). In addition, a desired target image pat- 
tern is provided, similar to the target pattern 106 of con- 
ventional MBOPC, as well as tolerance criteria, including a 
deviation tolerance, an image gradient tolerance and a 
curvature deviation tolerance. The edges of the primary 
mask features and the target image pattern are frag- 
mented into edge segments as in conventional OPC algo- 
rithms (Block 602). The image projected by the mask for 
the target lithographic process is simulated, including im- 
age intensity, and the image gradient (and other parame- 
ters as desired, for example, the curvature) at the selected 
(predetermined) comparison points (e.g. typically, but not 
limited to, one point at the center of each segment) (Block 
603). In order to compute the gradient and curvature, ad- 
ditional image values near the pre-determined evaluation 
point are typically required, for example, by computing 
values in a nine-point pattern around the evaluation 
point, which allows a finite difference approximation of 
the gradient and curvature, as described above (see Figure 
IE). 

[0052] | n accordance with the present invention, the deviation of 
the simulated image from the target image is estimated 
for regions between the predetermined evaluation points. 



If the estimated deviation exceeds a predetermined devia- 
tion tolerance, then a refinement of the relevant segment 
will be performed. The estimation of deviation is per- 
formed preferably by fitting a curve between adjacent 
evaluation points using the gradient and the curvature 
computed at the predetermined evaluation points. The 
gradient and curvature may be computed in a manner 
similar to that described previously, and are typically pro- 
vided by conventional OPC methods as previously dis- 
cussed. The computed gradient of the simulated image at 
the pre-selected points is first compared to a gradient 
tolerance (Block 605). If the gradient exceeds the gradient 
tolerance, then this indicates that a large deviation may 
exist between the target image and the simulated image, 
and therefore, at least one segment may need to be re- 
fined. The amount and nature of segment refinement is 
preferably determined by performing a curve fit (Block 
606), as described in more detail below. In the case that 
the computed gradient is zero (Block 613), it is possible 
that the image is at an optima (i.e. a maximum, a mini- 
mum, or a saddle point). In this case, the value of the cur- 
vature is used to determine whether a segment needs re- 
finement. If the curvature is also zero (Block 614) the seg- 



ment does not need refinement (Block 604) and the next 
segment is considered (Block 603). If the gradient is zero, 
and the curvature exceeds the curvature tolerance at the 
given predetermined evaluation point (Block 605), then 
that indicates that the segment may require refinement, 
and the method proceeds to estimation and evaluation of 
the deviation (Blocks 606-608). 
[0053] The process of segment refinement is shown in Blocks 
606 through 611. The first step in the process of refine- 
ment is to fit a curve using the image intensity, gradient 
and the curvature (the latter in case the gradient is zero) 
of adjacent evaluation points (Block 606). The next step is 
to compare the maximum deviation (i.e. distance) of the 
computed curve (i.e. the threshold intensity contour in the 
simulated image) from the target image (i.e. the threshold 
intensity shapes of the target image) for the correspond- 
ing segment. The curve fitting may be any method that 
fits a curve between two or more points, including, but 
not limited to, B-spline, polynomial approximation, circu- 
lar arc, cubic spline and Bezier curve. If the deviation of 
the curve exceeds the deviation (distance) tolerance (Block 
608), this indicates that the associated segment requires 
refinement, in accordance with the present invention. 



[0054] The corresponding segment for the mask and the target is 
then refined (Block 609) by breaking up each of the adja- 
cent segments into at least two segments. The newly cre- 
ated segments replace the original segments in the 
MBOPC algorithm. The image intensities are computed for 
these new segments (Block 610), in the same manner as in 
Block 603. If the image parameters (intensity, gradient 
and/or curvature) had been computed in a previous step 
and a segment has not gone through the refinement pro- 
cess, these values may be used in this step, without re- 
computing new values. 

[0055] This step is comparable to the step described in Block 104 
in Figure 1. The re-computed image intensities for the 
new segments are compared to the target image intensi- 
ties (according to the intensity tolerance) and the mask 
edge segments are modified accordingly (Block 611). This 
step is comparable to the step described in Block 107 in 
Figure 1. The evaluation of each mask and image segment 
is continued until no more segments require refinement 
(Block 604) and then the output mask layout is generated 
(Block 612). Figure 6 can be implemented on a digital 
computer as shown in Figure 12. 

[0056] Figure 12 shows a digital computer, whose main parts are 



the following: a central processing unit (CPU) 1201 , one 
or more input/output (I/O) devices 1202 (such as a key- 
board, a mouse, a compact disk (CD) drive, etc.), a con- 
troller 1203, a display device 1208, a storage device 1209 
capable of reading and/or writing computer readable 
code, and a memory 1206 - all of which are connected, 
for example, by a bus or a communications network 
(1205). The present invention can be implemented as a 
computer program product stored on a computer readable 
medium 1207, such as a tape or CD. The computer pro- 
gram product contains instructions to implement the 
method according to the present invention on a digital 
computer as shown in Figure 12. However, it can also be 
implemented in a multiplicity of such a digital computer 
where the items as depicted above may reside in close 
physical proximity or distributed over a large geographic 
region and connected by a communications network. 
[0057] a preferred embodiment of the adaptive refinement of the 
present invention can be better understood by reference 
to FIGs. 7-11. FIG. 7 illustrates a target image 400 having 
a (straight-edged) threshold boundary 410 including a 
primary feature having a non-critical portion 701 and a 
narrower portion 703 having a critical dimension (CD). 



The edges of the primary target feature 400 are seg- 
mented into segments having segment endpoints 710 
through 732. Superimposed on the primary target feature 
400 is the (curve-edged) boundary 413 of a simulated im- 
age 411 resulting from illuminating the primary feature 
400. Predetermined evaluation points are selected within 
each segment of the target boundary (e.g. predetermined 
evaluation points 802 through 836). 

[0058] The current invention works in the following manner. Im- 
age intensity, gradient and the curvature is computed at 
each of the pre-selected comparison points (Block 603). If 
the image parameters (intensity, gradient and/or curva- 
ture) had been computed in a previous step and a seg- 
ment has not gone through the refinement process, these 
values may be used in this step, without re-computation. 

[0059] Referring to FIG. 8, segment A, extending between points 
710 and 712 has a pre-selected evaluation point 802, and 
a segment adjacent to segment A, namely segment B, 
which extends between points 712 and 714, has a pre- 
selected evaluation point 803. The gradient at the pre- 
selected points 802, 803 is computed, for example, as 
discussed above, by using a finite difference approxima- 
tion to the gradient of equation (1) (e.g. see equation (2)), 



and curvature of equation (3). Since the values of gradient 
and curvature are typically computed in current state of 
the art methods, no additional computation cost is in- 
curred over conventional methods to obtain the gradient 
and curvature quantities for the present invention. 
[0060] The image gradient is compared to a gradient tolerance 
(Block 605). The gradient tolerance value is pre-specified 
along with the input mask. The gradient tolerance de- 
pends on the lithographic process condition and the cor- 
responding technology that is used to manufacture the 
wafer. If the gradient exceeds the gradient tolerance, then 
a curve fit to the intensity at the evaluation points on the 
target image is performed (Block 606). Note that in the 
case that the gradient is zero, the curvature should be 
checked (Block 614) to evaluate the case of an optimum or 
"saddle point" in the image intensity surface. A high cur- 
vature value, even when the gradient is zero, will indicate 
a large deviation of the simulated image from the target 
image. The curvature tolerance is pre-specified along with 
the input mask. The curvature tolerance depends on the 
lithographic process condition and the corresponding 
technology that is used to manufacture the wafer. For ex- 
ample, for current 193 nm lithographic processes, with 



numerical aperture about 0.75 and a standard (circular) il- 
lumination source, a gradient tolerance of about 6 and a 
curvature tolerance of about 10 may be appropriate, de- 
pending upon the desired tolerance in critical dimension. 
If both the gradient (and curvature values, if necessary) 
are less than the gradient tolerance (and curvature toler- 
ance, if necessary), then no further action is needed for 
that segment, and the next segment is evaluated (Block 
604). 

[0061] if the image gradient is greater than the gradient toler- 
ance, or if the curvature is non-zero and exceeds a curva- 
ture tolerance, in accordance with the present invention, 
the segment will be refined. The first step in a preferred 
embodiment of the present invention is to estimate the 
image between the evaluation point of the current seg- 
ment and an adjacent segment evaluation point based on 
the existing simulated image at those evaluation points. 
This estimated image is preferably determined by per- 
forming a curve fit using the intensity gradient and curva- 
ture (Block 606). At this step, in a preferred embodiment, 
a Binomial spline (B-spline) curve fit is performed using 
the information contained in the image intensities, image 
gradient and the image curvature. Other curve fitting 



techniques could also be used to implement the present 
invention. In this description of the invention the curve- 
fitting is used to denote any curve fitting method includ- 
ing but not exclusive of B-spline, polygonal approxima- 
tion, circular arc, cubic spline and Bezier curve. 

[0062] A n example of curve fitting is shown in Figure 8. In Figure 
8, a portion of the target image is shown having two seg- 
ments A (connecting points 710 and 712) and B 
(connecting points 712 and 714). In this example, the 
predetermined evaluation points 802 and 804 are at the 
center points of the two adjacent segments A and B. The 
arrows 852 and 854 show the directions of the image 
gradients at the evaluation points 802 and 804 respec- 
tively. By using the image intensity values and the gradi- 
ents at the evaluation points 802, 804, a B-spline curve fit 
840 is obtained using techniques known in the art. 

[0063] After the curve is fitted (Block 606), the maximum devia- 
tion distance 901 (D ) of the curve from the target im- 

max 

age is determined (Block 607). The maximum deviation 
distance 901 is preferably measured orthogonally from 
the target image. The orthogonal distance is a preferred 
deviation distance, because it is a very good norm for 45 
degree edges and for edges parallel to the coordinate 



axes. The maximum deviation distance 901 is depicted in 
Figure 9. The point on the target image where the maxi- 
mum deviation 901 occurs is denoted as P 

max 

[0064] a s depicted in Figure 6, D 901 is checked against a 

max 

predefined deviation tolerance (Block 608). If D is less 

max 

than the deviation tolerance, no further action is required 
for that segment, in this current step, in accordance with 
the inventive method, and the next segment is consid- 
ered. The deviation tolerance is defined by the user to ob- 
tain a desired accuracy within practical speed require- 
ments (i.e. a speed/accuracy tradeoff). A smaller tolerance 
value will result in increased refinements leading to more 
segments causing more evaluations and more iterations. 
On the other hand a larger tolerance value may not refine 
any segments making the resulting layout no better than 
that generated by the current art. A typical deviation tol- 
erance is about 5 -10% of the Critical Dimension (CD). For 
example, for a critical dimension of lOOnm , the deviation 
tolerance is preferably in the range of 5-10 nm. 
[0065] if d is greater than the predefined tolerance for dis- 

max 

tance, the inventive methodology then generates sub- 
segments (Block 609). The existing adjacent segments are 
divided so that the point of maximum deviation P be- 

max 



comes a new evaluation point 902 of a new segment A2. 
The new segment A2 is created by dividing the segment A 
in which the point of maximum deviation P is located, 

max 

into two sub-segments Al, A2, having new endpoint 
1002, which is preferably located so that the new evalua- 
tion point 902 (which is coincident with P ) is at the 

max 

center of a new segment A2. A new evaluation point 1102 
is also located within the other newly formed sub- 
segment Al, preferably at the center of Al. 
[0066] if t he point of maximum deviation P 902 were inside 

max 

the segment B with endpoints 712 and 714, then segment 
B needs to be divided, as illustrated in FIG. 10A. The seg- 
ment B is divided into two segments, viz., Bl and B2. In 
this case, the new segment endpoint 1002 is so chosen 
that the point 902 lies at the center of the newly created 
segment Bl, as in FIG. 10B. 
[0067] | n y e t some other cases, the point of maximum deviation 
P 902 may lie too close to or on the end point 712 di- 

max 

viding the segment A and B. Referring now to FIG. 10C, in 
this case, refining either segment A or B may create a 
segment with zero or very small length. Therefore, a new 
segment is created by taking portions of both segment A 
and B. In FIG. 10C, the newly created segment AB 1 has end 



points 1002 and 1003 and it satisfies the requirement of 
minimum length for a segment. The point of maximum 
deviation P 902 lies at the center of segment AB'. Con- 

max 

sequently both segments A' (with endpoints 710 and 
1003) and B' (with endpoints 714 and 1002) are created 
from segments A and B respectively. 

[0068] N 0 te that in the above described process, it was assumed 
that each segment is long enough to be refined. If a seg- 
ment is initially too small or becomes too small in the 
process of refinement, it may not be refined any further 
and an error message may be produced warning the user 
about this problem. 

[0069] The image intensity is computed at each of the new evalu- 
ation points 1102, 902 of the newly created segments 
(Block 610). If the threshold intensity of the image at the 
new evaluation point is not within the deviation distance 
tolerance limit relative to the target image threshold, the 
newly created segments are moved away or towards the 
target image for correction (Block 611). The process is 
then repeated with the new segments Al, A2. At the new 
evaluation points 1102, 902, new image intensity values, 
gradients and image curvatures are computed. Then the 
steps from Block 603 through 611 is repeated until the 



tolerances are met for all sub-segments of the original 
segment A being evaluated. For example, in Figure 11, the 
gradients 1152, 1154 are obtained for new evaluation 
points 1102, 902 respectively. The resulting fitted curve 
1111 has a maximum deviation distance D 1104, and a 

max 

P located at point 1106. If D at point 1104 exceeds 

max max 

the deviation tolerance, then the sub-segmentation con- 
tinues until the specified deviation tolerance is met. The 
resulting mask layout that is output at the end of the in- 
ventive method will have additional segments where 
needed to improve accuracy, but with only marginal in- 
crease in computational effort over conventional MBOPC 
methods. 

[0070] Although the present invention and its advantages have 

been described in detail, it should be understood that var- 
ious changes, substitutions and alterations can be made 
herein without departing from the spirit and scope of the 
invention as defined by the appended claims. 

[OO 71 ] WHAT IS CLAIMED IS: 



